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多 边 形 板 高 频 振动 的 能 量 虚 源 法 研究 


代 成 浩 , 钟 强 , 陈 海 波 
(中 国 科 学 院 材 料 力 学 行为 与 设计 重点 实验 室 ,中 国 科 学 技术 大 学 近代 力学 系 ,230026 合肥 ) 


摘 要 :本 研究 的 目的 是 将 室内 声学 中 的 虚 源 法 引入 高 频 横向 点 激励 作用 下 多 边 形 板 的 能 量 辐射 
传递 模型 中 ,以 预测 结构 弯曲 振动 的 能 量 响应 ,其 核心 是 如 何 构 建 多 边 形 板 的 虚 源 云 场 。 在 多 边 形 
板 的 能 量 辐射 传递 模型 框架 中 ,以 实 源 的 镜像 虚 源 蔡 代 结构 的 边界 ,以 此 来 表示 边界 处 能 量 的 连续 
反射 。 基 于 直接 场 的 能 量 密度 和 功率 流 强 度 核 函数 ,结构 内 部 任意 一 点 的 能 量 可 由 实 源 产生 的 能 
量 和 无 限 多 的 镜像 源 产 生 的 能 量 营 加 得 到 。 对 直角 三 角形 板 和 方形 板 求解 ,获得 了 它们 的 能 量 密 
度 与 功率 流 强 度 分 布 特征 。 数 值 算 例 结果 与 能 量 辐 射 传递 法 以 及 解析 解 对 比 ,验证 了 模型 的 准确 
性 。 最 后 讨论 了 阻尼 、 频 率 和 虚 源 的 阶 数 对 计算 结果 收敛 性 的 影响 。 
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Research on energy image source method for high frequency 
vibration of polygonal plate 
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Department of Modern Mechanics of University of Science and Technology of China ,230026 Hefei, China) 


Abstract : The purpose of this study is to introduce the image source method in room acoustics into the ra- 
diative energy transfer model of a polygonal plate under high frequency transverse point excitation to pre- 
dict its energy response of structural bending vibrations. The key treatment is how to construct the image 
source cloud field of the polygonal plate. In the framework of the radiative energy transfer model for the po- 
lygonal plate ,the mirror image of the real source is used to replace the boundary of the structure so as to 
represent the continuous reflection of energy at the boundary. Based on the energy density and power flow 
intensity kernel functions of the direct field , the energy at any point inside the structure can be obtained by 
superposition of the energy generated by the real source and the energy generated by an infinite number of 
image sources. By solving the right triangular plate and rectangular plate, we can obtain their distribution 
characteristics of energy density and power flow intensity. Numerical results are compared with those solved 
by the radiative energy transfer method and the analytical solutions to verify the accuracy of the proposed 


model. Finally, the effects of damping, frequency and order of image sources on the convergence of the 
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computation results are discussed. 
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大 量 的 单元 数目 和 对 不 确定 性 的 高 敏感 度 是 使 
MHA BN JE (finite element method, FEM ) LH .边界 元 
( boundary element method , BEM ) “等 确定 性 分 析 方 
法 求解 高 频 问 题 时 的 主要 限制 ,为 克服 这 些 限 制 ,者 
干 能 量 法 已 被 提出 并 在 高 频 分 析 领 域 成 为 了 主流 。 
其 中 ,最 常用 的 方法 是 统计 能 量 分 析 ( statistical ener- 
gy anlysis, SEA)  , 它 假设 能 量 场 为 扩散 场 ,但 只 
能 预测 系统 的 平均 能 量 响 应 。 为 此 ,Le Bot 在 SEA 
的 基础 上 ,提出 了 能 够 预测 结构 内 部 能 量 密 度 分 布 
的 能 量 辐射 传递 法 (radiative energy transfer method, 
RETM)“”。 该 方法 源 于 声学 中 的 射线 追踪 方 
渗 号 ,在 二 维和 三 维系 统 中 振动 场 的 形式 分 别 为 柱 
面 波 和 球面 波 ,结构 内 部 任意 一 点 的 能 量 由 实 源 辐 
射 的 直接 场 能 量 和 边界 虚 源 辐射 的 反射 场 能 量 琶 加 
得 到 ,在 高 阻尼 下 的 混 啊 场 仍 具 有 较 好 的 分 析 精 
HE^ 。 由 于 振动 传导 法 以 及 相应 的 能 量 有 限 元 法 
( energy finite element method , EFEM ) 假设 振动 场 由 
EET 288 f JL, AE, fe — He ux — HE RR EP, 
RETM 的 求解 精度 要 高 于 振动 传导 法 “”。 

虚 源 法 作为 一 种 预测 结构 高 频 振 动 响应 的 方法 
多 年 来 广 受 关注 “”。 最 早 , 虚 源 法 主要 被 应 用 于 
室内 声学 领域 ”| ,其 特点 是 以 实 源 在 域外 的 镜像 
虚 源 蔡 代 结 构 的 边界 ,基于 波 的 传播 规律 进行 计算 ， 
不 会 受到 常规 数值 分 析 中 单元 和 模 态 数 的 限制 。 
1905 年 ,Gunda 等 ”采用 虚 源 法 得 到 了 简 支 梁 和 和 矩 
形 板 的 振动 响应 。2009 年 ,Cuenca 等 ”应 用 虚 源 法 
得 到 了 任意 简 支 多边形 板 的 精确 格林 孔 数 ,有 效 的 预 
测 了 多 边 形 板 的 中 高 频 振 动 响应 。2012 4E, Cuenca 
等 同样 应 用 虚 源 法 得 到 任意 边界 条 件 下 多 边 形 板 
的 格林 函数 ,并 与 有 限 元 方法 结果 进行 了 对 比 。 

能 量 辐 射 传递 法 (RETM ) 是 通过 定义 在 边界 的 
虚 源 来 描述 能 量 的 反射 , 它 已 经 被 证 明 在 均匀 各 问 
同性 介质 的 同一 波 场 中 满足 射线 的 镜面 反射 规 
律 “” ,因此 ,RETM 在 理论 上 能 够 通过 虚 源 法 来 实 
现 。 虚 源 法 不 需要 求解 边界 虚 源 积分 方程 , 虚 源 的 振 
幅 为 实 源 振 幅 与 边界 反射 系数 的 乘积 ,只 需 对 所 有 虚 
源 的 影响 进行 全 加 , 即 可 获得 结构 内 部 的 能 量 响 应 。 

本 人 研究 通过 建立 多 边 形 板 的 域外 虚 源 云 场 , 采 
用 镜面 反射 来 描述 能 量 在 边界 处 的 传递 特性 ,将 虚 


源 法 推广 到 多 边 形 板 的 能 量 贺 射 传递 模型 中 。 数 值 
算 例 以 三 角形 板 和 和 矩形 板 为 例 ,获得 了 横向 点 激励 
下 结构 的 能 量 密度 与 功率 流 强度 。 最 后 讨论 了 不 同 
阻尼 和 不 同 虚 源 阶 数 下 ,基于 能 量 的 虚 源 法 与 能 量 
贺 射 传递 法 结果 的 误差 随 频 率 的 变化 。 


1 板 的 振动 模型 


1.1 板 中 的 波 传播 


槛 回 简 谐 点 激励 作用 下 的 多 边 形 板 模型 中 , 板 
的 振动 主要 为 模 癌 的 弯曲 振动 。 基 于 薄板 假设 的 多 
边 形 板 的 控制 方程 如 式 (1) 所 示 A 

DV'w+poh P = Fe"8(x-x)8(y-yy) (1) 


式 中 :w 2w(x,y) ARIS ETE; D = Eh /(12(1 - 
B. ) ) 为 弯曲 刚度 ;po 为 密度 ;h ON ERE SF. 是 频率 为 
co 简 谐 激励 的 振幅 ;6 为 Dirac delta KZ, Bix E 
控制 方程 的 行 波 解 为 w = Ae Ie nte 代入 
式 (1) 即 可 获得 板 弯 曲 振 动 的 频 散 方程 为 
Dk; ~pohw =0 (2) 
通过 式 (2) 获 得 弯曲 波 数 k, = (ooh AD)”, 弯 
曲 波 群 速度 c, 22 (Do / (ph) ) 


1.2 能 量 密度 与 功率 流 强 度 


根据 参考 文献 124], 板 内 任意 一 点 的 振动 由 其 
所 镜像 的 虚 源 云 场 的 所 有 虚 源 的 影响 倒 加 得 到 , 即 
任意 一 点 的 振动 可 由 式 (3) 表 示 为 


1 十 oo 十 oo 
W0(XY,Y3X0 yo) E " P 之 [Le(x,y;xo Yo) lak) x 
gf*y 


[w (x, y;0 001a 256 onem (3) 
式 中 ;A, A, 为 周期 性 虚 源 云 场 的 单元 场 尺寸 ;e 为 
描述 虚 源 云 场 单元 的 函数 ;w, 为 无 限 板 的 位 移 解 ; 
,=2mm/As ,k=2nm/A,。 结 合 特定 几何 形状 板 ， 
式 (3) 等 价 于 模 态 且 加 。 直 角 三 角形 板 和 和 矩形 板 的 
位 移 表 达 式 如 式 (4) 和 式 (5) 所 示 。 


4 X ws [sin(E,x,)sin(k,y,) 
wt . ) = | m” 0 n4 0 
enn pL Mgr EY E) 


sin( k, (L, m yo) ) sin( E, CL, 
D( (kn +h,) - kj) 


“ku Jin Ceo) sin( k, y) 


(4) 
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WC x yixy yg) = 
4 二 oo +o sin( k, xy ) sin( k, yy ) sin( k, x) sin( À, y) 
LL D (kn +k) - Kj) 


ym=1 n=1 


(5) 
因此 , c ALS h d 2 B BÉ t us E X CO) 
所 示 ” 。 


2 2 * 2 2 * 
w=- 4Re[D( 7-5 PA MAE 十 
4 OX OX 9y Oy 
2 2 x 2 2 c 
D eA su eee Es) 
Ox ðy oxðy ðxðy 
ðw Ow. 
"a a | D 


功率 流 强 度 与 能 量 密度 比值 为 群 速度 ,如 式 
(7) Biz c 


[=c,W (7) 


2C9 能 量 辐射 传递 模型 


2.1: 虚 源 云 场 


多 边 形 板 的 虚 源 云 场 坐标 主要 是 实 源 以 结构 边 
界 汐 镜面 获得 的 。 虚 源 法 是 以 波 传播 法 为 基础 的 分 
Mtn ds ,通过 确定 激励 点 与 接收 点 之 间 的 传播 路 径 ， 
预测 结构 的 振动 响应 。 

图 1 所 示 为 多 边 形 板 的 一 阶 虚 源 云 场 ,5。 为 作 
用 在 板 内 的 实 源 。51、5; 、S3 是 So 分 别 通过 边界 B, 
BB, 一 次 镜像 的 虚 源 。 


图 1 多 边 形 板 一 阶 虚 源 云 场 
Fig. 1 


First-order image source cloud 
field of polygonal plate 

图 2 为 多 边 形 板 二 阶 虚 源 云 场 , 图 中 S 为 经 过 
n 次 镜像 的 第 i 个 虚 源 。51 为 So 以 B, 为 镜面 的 一 
阶 虚 源 , 同 理 5; 和 5; 为 So 以 边界 B, 和 B, 为 镜像 
所 镜像 的 一 阶 虚 源 。 二 阶 虚 源 则 为 一 阶 镜像 虚 源 以 
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其 边界 为 镜面 获得 。 


2 S 
SE 
S 
图 2 多边 形 板 二 阶 虚 源 云 场 


Fig.2 Second-order image source cloud field of polygonal plate 
镜像 的 虚 源 坐标 可 由 式 (8) 获得, 即 

CS; -P): (P'-P) (P'-P) 

IP' -PI IP' -PI 
(8) 
其 中 :S ”为 点 S; 镜像 的 虚 源 ;P' 和 PP 分别 为 边界 
的 左右 两 端点 。 
能 量 辐射 传递 模型 

在 传统 的 二 维 结构 的 能 量 辐射 传递 模型 中 , 结 

构 任 意 一 点 的 能 量 为 实 源 辐射 的 直接 场 能 量 与 边界 

反射 场 能 量 的 车 加 ,通过 在 边界 处 划分 虚 源 来 描述 

能 量 的 反射 。 图 3 为 点 激励 作用 下 多 边 形 板 的 能 量 

辐射 传递 模型 示意 图 。 


St! = -S +2P +2 


2.2 


Mp; 
图 3 多 边 形 板 能 量 辐射 传递 模型 
Fig.3 Radiative energy transfer model for polygonal plate 


结构 的 能 量 密 度 和 功率 流 强 度 由 式 (9) 所 示 的 


积分 方程 获得 , 即 
W.(M) = |p(S)G,(S,M) ds + | a (P) G, CP, M) cosp dP 
1,(M) = [pC S) H,CS,M)4S + | o, (P) H, CP, M) cost; dP 


(9) 
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式 中 ;a 为 波 型 ;po 为 输入 功率 ;o 为 边界 虚 源 强度 。 
G, AH, 为 能 量 密度 和 功率 流 强 度 的 核 函数 "1 , 且 


- mr 
e 


G,(S,M) = 
2 TIC T sy 
n (10) 
H,(S,M) 23 —u(0) 
2arf oy 


式 中 :m = qo7 c, N BÉ EIER A, u (0) 2 3 e I8] 
Eo o 可 由 求解 边界 虚 源 功率 流 平衡 方程 获得 ,如 
A1) Bron, 

uS) 
Nu 


| [ Rp ( S) H, CS, P.) d 


[Rur PAL, „P,) cose dP,; |n», (11) 


将 虚 源 法 引入 到 RETM, 可 由 域外 的 虚 源 替代 
传统 RETM 定义 在 边界 的 虚 源 。 结 构 内 部 任意 一 点 
的 能 量 为 域外 虚 源 的 辐射 能 量 与 实 产 辆 射 能 量 和 县 
加 所 此 前 , 虚 源 法 预测 板 的 振动 是 从 位 移 的 角度 获 
得 振动 的 格林 函数 解 ” ,本 研究 从 能 量 的 角度 ,使 
用 虚 源 法 获得 结构 高 频 振动 的 能 量 啊 应 ,如 式 (12 ) 
所 不 。 


l 


WÈM) = P,[G,(S,,M) + Y, M RAG,CS;,M)] 
nz i=l 


(12) 
TSR: W 为 能 量 密度 ; P, 为 输入 功率 ;nm 为 虚 源 的 阶 
Dedh, 为 每 一 阶 虚 源 的 个 数 ;Rs 为 边界 能 量 反 冉 系 
数 呵 功率 流 强 度 如 式 (13 ) 所 示 。 
ICI) = P,[H, (S), M)u(0$) + 


Y Y RoHS, M) U( 0%) ] (13) 

虚 源 辐射 能 量 的 范围 是 有 限 的 ,如 图 4 所 示 ， 

一 个 虚 源 所 辐射 的 区 域 是 其 与 镜像 边界 两 端 过 线 中 

间 区 域 , 超 过 此 区 域 , 则 能 量 不 能 在 此 边界 反射 到 区 

域外 的 接收 点 。 当 虚 源 重合 时 , 虚 源 的 辐射 区 域 不 
同 , 因 此 只 计算 一 次 即 可 。 


S 十 +5 
B, 
S+ g| + 


图 4 虚 源 辐射 的 有 效 区 域 


Fig.4 The effective region of image sour ceradiation 
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3 数值 模拟 验证 与 讨论 


3.1 能 量 密度 及 功率 流 强度 


3.1.1 直角 三 角形 板 
以 直角 边 分 别 为 1 m x1 m 的 三 角形 铝板 为 例 ， 
厚度 为 1 mm, 四 周 简 文 ,其 材料 参数 如 表 1 所 示 。 
实 源 以 边界 镜像 的 前 五 阶 虚 源 云 如 图 5 所 示 , 在 
x, 70.25 m,y, 20.25 m 处 作用 五 =1N 的 横向 点 激 
励 ,频率 为 8 000 Hz. 
表 1 材料 参数 
Tab.1 Material parameters 
材料 ”弹性 模 量 E/GPa ” 泊 松 比 w —— WEHÉp/(kg: m^) 
Us 71 0.3 2 700 


yim 


x/m 
图 5 直角 三 角形 板 虚 源 云 场 


Fig.5 Image source cloud field of right triangle plate 


图 6 为 不 同 阻尼 下 ,由 式 (12) 得 到 的 能 量 虚 源 法 
结果 与 式 (4) 解 析 解 的 能 量 密度 (参考 值 10-“J/m ) 
的 对 比 。 从 图 中 可 以 看 出 ,两 种 阻尼 下 能 量 虚 源 法 
与 解析 解 吻 合 较 好 ,阻尼 越 大 能 量 衰 减 越 快 。 解 析 
解 中 徘 近 边 界 处 的 能 量 密 度 振 荡 较 大 ,这 是 由 于 反 
射 场 与 直接 场 的 干涉 引起 的 。 能 量 虚 源 法 忽略 了 波 
与 波 之 间 的 干涉 ,所 以 能 量 密度 的 分 布 较为 平滑 。 

图 7 和 图 8 分 别 为 2 种 阻尼 下 ,能量 虚 源 法 与 
解析 解 的 功率 流 强度 和 能 量 等 值 线 结 果 的 对 比 。 从 
图 中 可 以 看 出 ,阻尼 0.1 和 0.05 下 的 能 量 虚 源 法 与 
解析 解 的 能 量 流 方向 基本 一 致 ,能 量 等 值 线 的 分 布 
也 比较 吻合 。 板 的 振动 能 量 是 由 点 激励 引起 并 向 四 
周 流动 。 从 能 量 的 流 线 方向 和 等 值 线 的 分 布 可 以 看 
出 ,在 激励 点 附近 , 板 的 能 量 主要 受 直接 场 控 制 , 越 
靠近 边界 ,反射 场 的 影响 越 大 。 巾 于波 的 干涉 ,解析 
解 的 能 量 在 靠近 边界 的 位 置 ,振荡 明显 增 大 ,而 能 量 
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虚 源 法 的 结果 则 相对 平滑 一 些 。 


zi 
> 
m. 
co 
© 
= (b) 阻尼 比 为 0.05 
加 Klo 直角 三 角形 板 的 能 量 分 布 
Fig.6 Energy distribution of right triangle plate (b) 解析 解 
QN » 图 8 阻尼 比 w=0.05 时 直角 三 角形 板 的 功率 流 
CN Fig.8 Power flows of right triangle plate 
O when damping ratio ņ =0. 05 
N 3.1.2 方形 板 
之 0s 本 节 以 长 宽 高 分 别 为 1 m x 1 m x 1 mm 的 方形 
>< 铝板 为 例 ,四 周 简 支 ,其 材料 参数 如 表 1 所 示 。 在 板 
z 中 心 处 作用 振幅 F = IN ,频率 为 8 000 Hz 的 横向 点 
= 激励 。 研 究 不 同 阻尼 下 ,方形 板 能 量 密度 与 功率 流 
O , 强度 分 布 特征 ,并 与 解析 解 进行 对 比 验 证 引入 虚 源 
法 后 的 结构 能 量 辐射 传递 模型 的 准确 性 。 
(a) 能 量 虚 源 法 6 
T 
+ + 4 
" + + + + + 
+ + + + + + + 
下 
ES EocRoRooR +|+| + p xe ox we 
i 村 
a 
s ToLoRLokoc 
+ + + 
-4 f 
xim 
"— 图 9 方形 板 虚 源 云 场 
apo E Fig.9 Image source cloud field of square plate 
图 7 阻尼 比 m=0.1 时 直角 三 角形 板 的 功率 流 pen AEN QUA QE 
Fig.7 Power flows of right triangle plate 图 9 为 方形 板 前 五 阶 虚 源 云 场 。 图 10 为 采用 
when damping ratio 7 =0. 1 前 五 阶 虚 源 得 到 的 能 量 虚 源 法 结果 与 解析 解 的 能 量 
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密度 结果 的 对 比 。 从 图 中 可 以 看 出 ,依据 方形 板 的 
前 五 阶 虚 源 得 到 的 能 量 啊 应 与 解析 解 吻 合 较 好 。 当 
阻尼 为 0.1 时 ,能 量 由 激励 点 处 癌 四 周 快速 衰减 。 


(a) 阻尼 比 为 0.1 


ChinaXiv 合 作 期 刊 as 3o a 


而 阻尼 为 0. 05 时 ,能 量 的 衰减 速度 就 相对 绥 慢 
cir 
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(b) 阻尼 比 为 0.05 


图 10 方形 板 的 能 量 分 布 
Fig. 10 Energy distribution of square plate 


二 图 11 和 图 12 为 方形 板 的 功率 流 强度 分 布 。 从 
图 中 可 以 看 出 ,能 量 虚 源 法 与 解析 解 的 能 量 流 线 的 
和 等 势 线 基本 一 致 。 与 直角 三 角形 板 的 能 量 响 
应 贰 果 有 着 同样 的 现象 ,由 于 忽略 了 波 和 波 之 间 的 
析 合 干涉 ,因此 能 量 虚 源 法 的 结果 较为 平滑 。 而 解 


202211 


chinaXiv 


(a) 能 量 虚 源 法 


图 11 阻尼 比 7 =0.1 时 方形 板 的 功率 流 强度 


Power flows of square plate when damping ratio 7 =0. 1 


Fig. 11 
3.2 计算 收敛 性 分 析 


本 市 以 上 节 方 形 板 模型 为 例 ,讨论 不 同 阻尼 下 
的 频率 和 虚 源 阶 数 对 局 部 能 量 响应 与 平均 能 量 啊 应 
收敛 性 的 影响 。 图 13 为 (0.8 m,0. 8 m) 和 (0.4m， 
0.4m) 处 的 局 部 能 量 啊 应 随 频率 的 变化 。 从 图 中 可 
以 看 出 , 当 激 励 点 取 在 板 中 间 时 , 随 厦 频率 的 增加 ， 
两 种 阻尼 下 的 能 量 虚 源 法 结果 与 RETM 解 基本 吻 
合 , 且 均 逐 渐 收 敛 于 解析 解 。 由 于 (0.8 m,0.8 m) 点 
处 更 接近 边界 ,其 解析 解 的 振荡 幅度 相对 大 些 。 随 
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析 结 果 考 虑 了 直接 场 和 反射 场 之 间 的 干涉 效应 ， 
此 能 量 在 边界 处 的 振荡 较为 明显 。 由 于 方形 板 的 结 
构 具 有 对 称 性 ,激励 点 取 在 板 的 中 心 , 因 此 能 量 密度 
的 结果 也 呈现 出 对 称 现象 。 


x/m 


(b) 解析 解 


着 频率 的 增加 ,由 激励 位 置 向 边界 的 能 量 衰减 程度 
逐渐 增 大 ,阻尼 越 大 ,此 衰减 幅度 也 越 大 。 

为 体现 虚 源 的 数量 对 结果 收敛 性 的 影响 ,以 
RETM 的 结果 为 参考 值 ,分 析 能 量 虚 源 法 与 RETM 
的 能 量 结 果 的 偏差 随 频 率 和 虚 源 阶 数 的 变化 。 图 
14 为 采用 能 量 虚 源 法 与 传统 的 RETM 计算 的 
(0.4 m,0. 4 m) 处 能 量 响 应 的 偏差 分 析 结 果 。 从 图 
中 可 以 看 出 随 着 频率 的 增加 ,能 量 虚 源 法 逐渐 收敛 
于 RETM 解 。 当 参与 计算 的 虚 源 阶 数 越 多 时 ,收敛 
越 快 。 当 虚 源 阶 数 一 定时 ,阻尼 越 大 ,结果 随 频率 的 
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M 
增加 ,收敛 速度 越 快 。 当 频率 一 定 ,阻尼 越 大 时 ,收敛 ”结构 精确 的 能 量 响应 所 需要 的 虚 源 数 数量 相对 
于 RETM 的 结果 值 所 需要 参与 计算 的 虚 源 数 越 少 。 较 少 。 

此 , 当 激 励 频率 和 阻尼 均 较 大 时 ,采用 能 量 虚 源 法 获得 


1.0 o 


x/m 
Y (a) 能 量 虚 源 法 (b) 解析 解 
a 图 12 BEEE y =0. 05 时 方形 板 的 功率 流 
CO Fig. 12 Power flows of square plate when damping ratio ņ = 0. 05 
C 110 110 
© (0.4 m,0.4 m)RETM (0.4 m,0.4 m)RETM 
C 100 (0.4 m,0.4 m) 能 量 虚 源 法 100 (0.4 m,0.4 m) 能 量 虚 源 法 
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© 10 10° 10° 10 10' 10 10 10° 
= f/Hz f /Hz 
= (a) 阻尼 比 为 0.1 (b) 阻尼 比 为 0.05 
O 图 13 随 频率 变化 的 局 部 能 量 响应 


Fig. 13 Local energy response that varies with frequency 
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(a) 阻尼 比 为 0.01 (b) 阻尼 比 为 0.03 


图 14 〈0.4m,0.4m) 处 能 量 虚 源 法 与 RETM 局 部 能 量 啊 应 偏差 分 析 


Fig. 14 Analysis of local energy response deviation between image source method and RETM at (0.4 m,0.4 m) 


图 15 为 采用 能 量 虚 源 法 与 RETM 计算 得 到 的 “平均 能 量 啊 应 的 误差 随 频率 变化 的 趋势 。 从 图 中 可 
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以 看 出 平均 能 量 啊 应 变化 的 趋势 与 图 14 中 的 局 部 

能 量 啊 应 偏差 趋势 基本 一 致 。 当 频率 一 定时 ,参与 

计算 的 虚 源 数量 越 多 , 结 末 越 精确 , 且 当 虚 源 数量 增 
100 
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(a) 阻尼 比 为 0.01 
图 15 平均 
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学 学 Od = 了 第 39 卷 
加 到 一 定 程度 时 ,结果 逐渐 收敛 。 当 能 量 虚 源 法 的 
结果 收敛 时 ,与 RETM 的 平均 能 量 响 应 还 是 存在 略 


微 的 偏差 。 
100 
一 1 阶 e 5 阶 
-—— 2 阶 9 6 阶 
80 ---3 阶 9 7 阶 
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(b) 阻尼 比 为 0.03 
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Fig. 15. Average energy response deviation analysis 


OO 由 于 能 量 虚 源 法 采用 的 是 镜面 反射 模型 ,而 传统 
的 RETM 边界 处 采用 漫 反 射 模型 描述 边界 处 能 量 的 
反对 ,为 体现 反射 模型 的 影响 ,在 方形 板 的 (0.25 m, 
0.25 m) 处 作用 一 幅 值 为 1 的 简 谐 点 激励 ,结构 阻尼 
79:1, 

^U 16 为 当 阻尼 比 为 0.1 时 ,在 8 000 Hz 和 


1.0140 "40 


ym 


x/m 
(a) 频率 为 8 000 Hz 


图 16 


12 000 Hz 的 频率 下 采用 RETM .能量 虚 源 法 和 解析 
方法 获得 的 方形 板 的 能 量 等 势 线 结果 对 比 。 从 图 中 
可 以 看 出 在 不 同 的 频率 下 ,采用 能 量 虚 源 法 得 到 的 
能 量 分 布 结果 相 较 于 RETM 与 解析 解 更 加 吻合 一 
些 。 这 是 由 于 采用 镜面 反射 比 漫 反 射 更 符合 真实 波 
场 的 反射 规律 。 

7 30 一 一 T 

| um m c feur 


ud  —]—— g 


— — REIM 


x/m 
(b) 频率 为 12 000 Hz 


能 量 等 势 线 


Fig. 16 Energy contours 


eh 论 


本 研究 以 多 边 形 板 为 研究 对 象 , 建 立 了 通过 虚 
源 法 实现 的 多 边 形 板 能 量 辐射 传递 模型 。 以 直角 三 
角形 板 和 方形 板 为 例 , 建 立 了 其 虚 源 云 场 ,获得 模型 
在 横向 点 激励 作用 下 的 弯曲 波 能 量 密度 和 功率 流 强 
度 分 布 特征 。 在 不 同 的 阻尼 下 ,数值 结果 与 解析 解 


和 传统 的 RETM 结果 进行 了 对 比 ,验证 了 模型 的 准 
确 性 并 分 析 了 频率 对 局 部 能 量 响应 的 影响 。 最 后 以 
方形 板 为 例 ,讨论 了 能 量 虚 源 法 与 RETM 的 偏差 随 
频率 .阻尼 及 虚 源 的 阶 数 的 变化 ,所 得 结论 如 下 。 

1) 直角 三 角形 板 和 方形 板 的 能 量 密度 与 解析 解 
在 不 同 阻尼 情况 下 均 吻 合 较 好 ,功率 流 强度 分 布 特 
征 也 与 解析 解 较为 一 致 。 

2) 能 量 虚 源 法 的 结果 与 解析 解 随 着 频率 的 增加 
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吻合 程度 越 好 。 随 着 频率 的 增加 ,其 计算 结果 逐渐 
收敛 于 RETM 解 , 当 参与 计算 的 虚 源 阶 数 越 多 时 , 收 
敛 越 快 。 阻 尼 越 大 时 ,收敛 于 RETM 解 所 需要 参与 
计算 的 虚 源 数 越 少 。 

3 ) 能 量 虚 源 法 相 比 于 传统 的 RETM 与 解析 解 更 
为 吻合 ,但 与 RETM 的 偏差 并 不 大 , 属于 可 接受 
范围 。 
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